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I review what we know about the "sphaleron rate" , which is the efficiency of 

O ■ baryon number violation at high temperatures T ~ 100 GeV in the Standard 

Q^ ' Model. The leading behavior at weak coupling in the symmetric phase is known 

Ft ' accurately; V = (10.7 ± .7){g'^T^ /m\,)\og{mu / g'^T)al,T'^ . At realistic values of 

the coupling our accuracy is worse. We also now have the tools to determine the 
Vi , rate nonperturbatively in the broken electroweak phase; the sphaleron rate there 

Cm ' is slower than perturbative estimates. 

T-H ■ 
K^ \ 1 Introduction 

\^ • The Universe is filled with matter, and virtually no antimatter. This "unusual" 

^T I situation can only be explained without appealing to the initial conditions of 

the Universe if, at some early epoch in the history of the universe, baryon 
number was not a conserved quantity, 
^s . As a matteiL of fact, baryon number is not a conserved quantity in the 

standard model El. Furthermore, while its violation under ordinary conditions 
is pitifully inefficient, that ceases to be true at very high temperatures, of order 
the weak scale, where electroweak symmetry is restored. These facts are the 
j^i backdrop for the subject of electroweak baryogenesis, which attempts to use 

them to explain why the baryon number density of the current universe is what 
it is. 

In this talk I will not discuss baryogenesis. Rather my emphasis is on 
strengthening its foundations by investigating more accurately exactly how 
^ I efficiently baryon number is violated under hot conditions. Besides the obvious 

application to the study of baryogenesis, this is also a useful thing to do because 
it forces us to develop tools for dealing with the infrared physics of hot Yang- 
Mills theory. Also it is the only part of a baryogenesis calculation which is 
generic to all extensions of the standard model, because it only depends on 
the gauge group, and only very weakly on the Higgs sector. In fact, for much 
of the talk I will neglect the Higgs fields altogether and just study Yang-Mills 
theory. 

1.1 What we want to know 

Before going any further I should fix notation and explain what I will measure. 
The anomaly equation relates baryon number to the Chern-Simons number of 
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the SU(2) weak fields, through 

iN^it) ^Ncs^^J dt' j cPxE^Btix, t') , (1) 

where E and B are the SU(2) electric and magnetic fields, and I normalize so 
the gauge field A has units of inverse length and g^ appears in the denomi- 
nator in the Lagrangian. The constant of integration from the indefinite time 
integral is fixed by requiring A^'cs to be an integer for a vacuum configuration. 
iVcs is called the Chern-Simons number. Because magnetic fields are always 
transverse (Gauss' Law for magnetism), the evolution of Chern-Simons number 
depends on the physics of the transverse sector. 

Having defined Chern-Simons number I can define its diffusion constant, 

r.hm hmii^M^l^JYcsM!), (2) 

where the angular brackets () mean an average is taken over the thermal density 
matrix. F is often referred to as the "sphaleron rate." The reason we care 
about it is that there is a fluctuation dissipation relation between it and the 
relaxation rate for a chemical potential for baryon number. I will not discuss 
this in detail, see insteadS'crQ. I also comment that for iVcs to diffuse requires 
nonperturbative physics, and nonperturbative physics is only unsunpressed, at 
high temperatures and weak coupling, on length scales > (l/g^T) H. 

It is interesting to know T in two regimes. The first is the electroweak 
symmetric phase. Almost all baryogenesis mechanisms will give a final baryon 
number directly proportional to its value here. For this reason we would like 
to know it with some accuracy here, which makes the calculation tricky. The 
other regime where we want to know F is in the broken electroweak phase, 
immediately after the electroweak phase transition, which is to say, right after 
the baryons were allegedly produced. In this case what we want to know is, are 
the baryons safe, or will they subsequently be destroyed? Here the strength 
of the phase transition is important; the real question is, how strong must the 
phase transition be to prevent the baryons from getting destroyed? To answer 
this question with descent resolution we only need to know F to within ±1 in 
the exponent; however F is exponentially small and perturbation theory is not 
yet reliable, which will make this calculation tricky as well. 

1.2 Approxim,ations I will need 

Determining F requires determining unequal Minkowski time correlators at 
finite temperature in a quantum field theory. Furthermore, if the answer is to 



be nontrivial the field theory must be showing nonperturbative physics. No 
one knows how to do this directly. Therefore I will obviously need to make 
some approximations. 

What saves us is that the SU(2) sector is weakly coupled. This allows 
two key approximations which make the problem tractable. First, the infrared 
behavior of the theory is classical up to parametrically suppressed corrections. 
Second, the ultraviolet behavior of the theory is perturbative. Here, infrared 
means k <C ttT, while ultraviolet means k ^ a^jT. When the coupling really 
is weak, there is an overlap between these two regimes, and every degree of 
freedom can be treated with one approximation or the other. Then, one can 
integrate out those degrees of freedom which are perturbative, and treat the 
remaining, classical theory nonperturbatively on the lattice. 

In what follows I will first discuss what we learn by treating perturbatively 
and integrating out everything we can. Then I will step back and only integrate 
out the highest k modes analytically, leaving a larger and more inclusive theory 
for numerical work. Finally I discuss what we can do in the broken phase. 
The complete details for these three approaches can be found in Q, 0, and H 
respectively. 

2 Leading log 

Dietrich Bodeker has shown that it is possible to integrate out all degrees of 
freedom with momentum scale k > g^T\og{l/g), and that doing so produces 
an effective theory for the remaining k '^ g^T degrees of freedom which is 
classical Yang-Mills theory under Langevin dynamics B. The xihysical origin 
of this effective theory has been discussed by Arnold's group tj. Integrating 
out thc_inodes with k > gT gives the well known hard thermal loop effective 
IheoryEJ. The behavior of the infrared modes in this theory is overdamped 
E3, which just follows from Lenz's law and the fact that the plasma is highly 
conducting. The conductivity is k dependent on scales shorter than some mean 
collision length, which in an abelian theory is the large angle scattering length. 
However in a nonabelian theory a particle's charge is changed by scattering. 
The mean length for a particle to travel before its charge is randomized is 
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SO on scales longer than this the strength of damping is k independent. There- 
fore, in the approximation that the scale 1/g^T is well separated from the 
scale 1 / {g^T \og{l / g)) , the infrared fields obey Langevin dynamics on long 
time scales. 



lattice spacing a 


Volume 


Langevin time 


k'± statistical error 


2/ig^T 


{^/g'Tf 


290000a^ 


10.44 ±0.23 


2/ig^T 


[m/g'Tf 


49500a^ 


10.30 ±0.21 


2/bg'T 


{m/g'Tf 


21000a^ 


10.70 ±0.67 


2/7g'T 


{IQ/g'Tf 


42000a^ 


10.25 ±0.79 



Table 1: Results for k' at three lattice spacings and two lattice volumes. The results show 
excellent spacing and volume independence. 



Langevin dynamics have two nice features. First, the Langevin dynamics 
of 3-D Yang-Mills theory are free of UV problems, and a zero lattice spacing 
limit exists. Second, Langevin dynamics are very easy to put on the lattice. 
Therefore the emphasis should be on controlling systematics, such as 

1. the thermodynamic match between lattice and continuum, 

2. the match between lattice and continuum Langevin time scales, 

3. topological definition of iVcs, and 

4. the large volume and long time limits. 

All of_these systematics can be controlled. The first is discussed intj, the 
second inQ, and the third inB. A volume S/g^T on a side is large enough to 
achieve the large volume limit, so I use a volume IQ/g'^T on a side for overkill. 
The failure to achieve the large time limit is refiected in the statistical error 
bars. 

The results, which show beautiful lattice spacing independence, are pre- 
sented in Table |l|, which gives the coefficient k! for the equation 
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These results settle the question, "What is the sphaleron rate in the symmetric 
phase, in the extreme weak coupling limit?" 

3 Beyond the leading log 

In the last section I determined the coefficient of the leading log behavior with 
very good precision. The problem is the (+0(1)) appearing in Eq. (0). How 
well does the expansion in \og(\/ g) converge? 

The answer is, probably very poorly. To see this, compare the free path 
for color randomization, Zmax = 2iT/g^T\og{l/g), to the size of a box for which 




Figure 1: Sphaleron rate plotted against inverse HTL strength, for the lattice theory plus 
"particles" used to generate the hard thermal loops. The dashed fit is a straight line, while 
the dotted fit incorporates the known logarithmic dependence on mn found in the previous 
section. 



the large volume limit has already been reached, ^/ g^T. There is not a large 
separation between these scales. In fact, it is not clear whether Zmax is smaller 
than the scale characterizing nonperturbative physics, which must after all be 
well shorter than the dimension of a box which shows large volume behavior. 
The problem is that Bodeker's approach requires integrating out modes for 
which the perturbative treatment may not be very reliable. To test this, and 
to try to determine the sphaleron rate beyond leading log, we need to integrate 
out less, and make the numerical model include the gT as well as g^T scales. 
An effective action for the theory with the k ^ T modes integrated out 
is known Ej, and goes by the name of the hard thermal loop (HTL) effective 
action. It is nonlocal, which is not surprising, since its construction involves 
integrating out propagating degrees of freedom in a Minkowski theory. Unfor- 
tunately nonlocality is very problematic for a numerical implementation. 
i— I A solution to this problem was proposed a few years ago by Hu and Miiller 
lla. The idea is that, rather than add the HTL action itself, one adds a set of 
local degrees of freedom which, if integrated out, would also yield the HTL ef- 



fective action. Since the HTL action represents the propagation of a set of high 
momentum, charged particles, what we add are a bunch of high momentum, 
charged classical particles. For the idea to work it is necessary to add particles 
in such a way that the numerical model retains an exact gauge invariance, and 
possesses a conserved energy and phase space measure so that thermodynamic 
averages are well defined. We present an implementation which satisfies these 
conditions inQ, and refer the reader there for the (quite complicated) details. 
With a model which reproduces the HTL resummed IR physics in hand, 
it is possible to test Arnold, Son, and Yaffe's claim that the IR dynamics are 
overdamped, directly. The results are shown in Figure |^, which shows that 
r does vanish linearly as mj~, is increased. The data are not good enough 
to show unambiguously whether or not Bodeker's log is present. Fitting the 
data assuming it is, the (+0(1)) in Eq. (||) turns out to be about 3.6, which 
indicates that the expansion in log{mu/g^T) is not a very good one. For 
the physical value of mj^ = {ll/6)g^T^ and g^ ^ 0.4, the sphaleron rate is 
(20 — 25)a^r'*, with systematics dominated errors of the order of 30%. The 
main remaining problems to be addressed involve lattice spacing effects and 
the interactions of the added "particle" degrees of freedom with the most UV 
lattice modes. 



4 The broken electroweak phase 

In the previous two sections the approach has been to find a numerical system 
which has the same physics as thermal Yang-Mills theory, and then to evolve 
it and measure directly the correlator, Eq. (0), which tells how efficiently 
baryon number is violated. However, this approach fails completely in the 
broken electroweak phase, because the rate of topological transitions is so small 
that no reasonable amount of numerical evolution would see any transitions 
at all. Another alternative, perturbation theory, is not very reliable close to 
the electroweak phase transition. We know for instance that the one loop and 
two loop effective potentials give quite different answers for the strength of the 
transition, and no one knows how to compute the sphaleron rate beyond the 
one loop level. Some other technique, nonperturbative but not strictly real 
time, is needed. 

The reason for the suppression of the sphaleron rate in the broken phase 
is shown, in cartoon form, in Figure 0. There is a free energy barrier between 
minima, meaning that almost none of the weight of the thermal ensemble lies 
in states intermediate between vacua. The figure also suggests how I will 
determine the sphaleron rate in the broken phase. I should define A'^cs (or 
some appropriate observable) on the lattice, and measure how the free energy 
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Figure 2: "Cartoon" of the free energy dependence on A^cs- 



depends on it. This is not enough to give the real time rate, but with some 
more work one can turn the height of the barrier into the real time rate. 

4-1 Defining Ncs 

I begin by defining Chern-Simons number on the lattice. The obvious approach 
is to use the same definition as in the continuum, Eq. (|l|). Note that the 
integral over "time" in that equation could really be an integral along any 
path through the space of configurations, not just one generated by Hamilton's 
equations. In particular one can fix the constant of integration by having the 
path begin or end at a vacuum configuration. 

There is a problem on the lattice, which is that no lattice implementation 
of EfBf is exactly a total derivative. Therefore iVcs defined through Eq. 
(|l|) and implemented on the lattice would depend on the path chosen. The 
resolution is to choose a unique and particularly sensible path, the gradient 
flow (cooling) path, that is, the path through configuration space along which 
the energy falls most rapidly. In continuum notation the path (parameterized 
by a cooling time r) is given by 

dA{x,T) _ dH 

dT ~ dA{x) ' ^ ' 

where H is the Hamiltonian and all indices have been suppressed. Besides 
being unique, this path also has the benefit that it moves quickly towards 
configurations where the gauge fields are smooth. This minimizes the impact of 
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lattice artifacts, which were for instance responsible for EfBf not being a total 
derivative. Further, the path automatically goes to a vacuum configuration. 

This definition of iVcs has two problems, both easily resolved. First, TVcs 
is UV poorly behaved. For instance, its mean squared value diverges as V/a, 
with V the physical volume and a the lattice spacing (or other regulator) . This 
is resolved by using not A^cs but the Chern-Simons number of a configuration 
after an initial length tq of gradient flow. Our measurable is then dependent 
on an unphysical parameter rp, but it is UV finite, and physical measurables 
such as r will be tq independent in the end. 

The second problem is that performing gradient flow down to the vacuum 
is intensely numerically expensive; yet we will need to do so thousands of 
times to determine the free energy distribution. This problem is solved by 
blocking. Gradient flow destroys information, and in particular it destroys 
almost all the UV information; so nothing is lost by blocking after some modest 
amount of gradient flow. The numerical savings are immense, and (if we use 
an 0{a?) improved lattice Hamiltonian and implementation of EfBf) almost 
no accuracy is lost. 

Using this definition of Nq^^ it is possible to determine the free energy 
(probability distribution) as a function of Nqc^ by standard multicanonical 
Monte-Carlo techniques. A sample result is shown in Figure 0. 

4. 2 Turning probabilities into rates 

One cannot read off the sphaleron rate from Figure alone; in fact the height 
of the barrier in the figure depends on an arbitrary parameter tq. To get F 
from the figure we need to know 



{N) 



dNcs 



dt 



), (6) 

Afcs=0.5/ 



the mean rate at which Ncs is changing during a crossing of the barrier. 
Multiplying the probability density at the top of the barrier by (N) turns the 
probability density into a probability fiux per unit time. 

It is straightforward to measure (N) numerically. First, we use multi- 
canonical means to get a sample of configurations with Nqq ~ 0.5. Then, for 
each we draw momenta randomly from the thermal ensemble and perform a 
very short period of Hamiltonian evolution, measuring A'^cs before and after. 
Then \dNcs/dt\ is approximated by |iVcs(0) - Ncsi^t)] / (St) , and {N) is the 
average of this over the sample. Also, one must divide by the volume used in 
the lattice simulation, to convert the rate of topological transitions to the rate 
per unit volume. 
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Figure 3: Free energy as a function of A^cSi ^i the critical temperature in a {16 /gT) 
cubic box and using tq = 3.6/(g^T)^, when the ratio of the Higgs sclf-coupUng to the gauge 
coupling was \/g^ = 0.039. The upper curve is the broken phase value and the lower curve 
is the symmetric phase value. 



However, F does not equal the probability flux per unit time over the 
barrier; it is how often one goes from being in one topological vacuum to being 
in another. It is possible to cross the barrier several times on the way from one 
minimum to its neighbor, or to cross an even number of times and return to the 
starting vacuum. This leads to a correction called the "dynamical prefactor," 
which is the ratio of true topological vacuum changes to crossings of the top 
of the barrier. To compute it, we use multicanonical means to get a sample of 
A^cs = 0-5 configurations. Then each is evolved under Hamiltonian dynamics, 
both forward and backwards in time, until it settles in a topological vacuum. 
The dynamical prefactor is 



Prefactor = 



V '- 

^-^ # crossings 

sample 



(ATVcs)' , 



(7) 



where ATVcs is the difference in iVcs between the starting and ending vacua. 
It is ±1 if there were and odd number of iVcs ~ 0.5 crossings and if there 
were an even number; we never observe prompt crossings from one topological 




Figure 4: Chern-Simons number diffusion constant T immediately after the electroweak 
phase transition as a function of \/g^, which at tree level equals m'jj/8m?y. The solid line 
is a perturbative estimate, the line at the top of the figure is the symmetric phase rate. 



minimum to another which is not its immediate neighbor. 

The hard thermal loops appear in the dynamical prefactor, which is para- 
metrically of order {g T^/mj~,)log{mi:)/g^T). However, using the techniques 
of the last section to include the HTL effects in the calculation of the prefactor 
shows that, for realistic values of mfy, the importance of HTL's is weak. This 
is expected, or at least we should expect that the dependence is weaker than 
in the symmetric phase, because broken phase baryon number violation should 
be mediated by a spatially smaller configuration, involving higher frequency 
modes which are less overdamped. 

The final result for F in the broken electroweak phase, at the electroweak 
phase transition temperature and for a range of scalar self-couplings, is plotted 
in Figure ^ which also compares it to a perturbative result based on a tazp loop 
potential and the zero mode calculations from Carson and McLerran E3. The 
actual rate is substantially but not drastically slower than the perturbative 
estimate. For comparison, the value needed to avoid baryon number washout 
after the transition, in the standard cosmology, is F ~ lO^^a^T^. 
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5 Conclusion 

Tools now exist to calculate the baryon number violation rate in both the 
symmetric and broken electroweak phases. In the symmetric phase the rate 
behaves parametrically as a^, with a logarithmic correction found by Bodeker 
which is numerically small. The rate at a realistic m|, is 20 — 25a^T'*, with 
systematic errors, estimated to be of order 30%, dominating statistical errors. 
In the broken phase the rate is smaller than a perturbative estimate, but still 
too large to save baryogenesis in the minimal standard model. 
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